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NiO is a good candidate for ultrafast magnetic switching because of its large spin density, antiferromagnetic 
order, and clearly separated intragap states. In order to detect and monitor the switching dynamics, we develop 
a systematic approach to study optical second harmonic generation (SHG) in NiO, both at the (001) surface 
and in the bulk. In our calculations NiO is modeled as a doubly embedded cluster. All intragap d-states of the 
bulk and the (001) surface are obtained with highly-correlational quantum chemistry and propagated in time 
under the influence of a static magnetic field and a laser pulse. We find that demagnetization and switching can 
be best achieved in a subpicosecond regime with linearly rather than circularly polarized light. We also show 
the importance of including an external magnetic field in order to distinguish spin-up and spin-down states and 
the necessity of including magnetic-dipole transitions in order to realize the A-process in the centrosymmetric 
bulk. Having already shown the effects of phonons in the SHG for the bulk NiO within the frozen-phonon 
approximation, and following the same trail of thoughts, we discuss the role of phonons in a fully quantized 
picture as a symmetry-lowering mechanism in the switching scenario and investigate the electronic and lattice 
temperature effects. 

PACS numbers: 78.68.+m, 78.47.+p, 73.20.-r, 78.20.Ls 



I. INTRODUCTION 

In recent decades there has been a continuous strain to min- 
imize as much as possible the time needed to record data on 
magnetic materials. At the same time, the tremendous in- 
crease in storage density and read-write speed in magnetic 
storage media is reaching its physical limits. Ever since the 
light-induced demagnetization of ferromagnets was discov- 
ered', several light-driven scenarios and mechanisms have 
been proposed^'^, and it has already been demonstrated that 
exploiting the ultrafast electron-photon interaction can lead to 
subpicosend dynamics''. Furthermore, it has been experimen- 
tally assessed that the electronic excitations can non-thermally 
control spin dynamics^ and magnetic phases^ in magnetic 
materials. Laser manipulation has been well understood in 
atomic and molecular systems, which stimulates the demand 
for an implementation of these mechanisms in practical de- 
vice applications. Thus, the investigation of electron and spin 
dynamics in the bulk, surfaces, and nanostructures of nonmag- 
netic as well as magnetic materials has attracted considerable 
interest in recent years. 

Nonlinear optics is more selective than nonlocal linear op- 
tics, making it an ideal tool for investigating antiferromagnets. 
NiO is a good candidate for such a scenario due to its large 
spin density, antiferromagnetic (AF) order, and clearly sepa- 
rated intragap states. Although bulk NiO possesses a center 
of inversion, it generates a second harmonic signal^. In the 
literature there are todate four main explanations for its origin 
(i) spin-orbit coupling (SOC), (ii) a signal that mainly results 
from the surface, (iii) inclusion of higher order transitions, 
and (iv) lowering of the crystallographic symmetry due to lo- 
cal distortions or phonons. Recent works^-^ from our group 
showed that explanations (i) and (iii), especially the latter, are 
the most probable candidates for the physics behind. The first 
two explanations are also demonstrated by a recent experi- 



mentiS. 

In this paper, we review our studies on electronic theory 
of ultrafast spin dynamics in recent years. The remainder 
of this paper is organized as follows. In Sec. 2, we intro- 
duce quantum chemistry methods used in the calculations. In 
Sec. 3, we include SOC. In Sec. 4, second harmonic gen- 
eration (SHG) is calculated. Symmetry analysis of the con- 
tributions to the second order susceptibility tensors (x'-^'^') is 
performed. In Sec. 5, we show the effects of phonons as a 
symmetry-lowering mechanism in the switching scenario. In 
Sec. 6, we present our final results for realization of spin ma- 
nipulation. Finally, in Sec. 7, we present our summary and 
outlook for future work. 



II. QUANTUM CHEMISTRY 

One can use either the real or the momentum space ap- 
proach to model both the surface and the bulk of NiO. Their 
major difference is that only the first one is able to find the 
localized intragap c/-states of the Ni ion. These dispersion- 
less states are experimentally confirmed for the surface iij^^ 
and for the bulk'^. Although some previous calculations 
with extended local density approximation (LDA-H-) within 
the Hubbard-I approach^ and within the GW approxima- 
jjQj^i5ji6ji7 gjyg better results for the gap than the real space 
approach, and a recent calculation combining LDA and dy- 
namical mean-field theory provides an accurate band structure 
of NiOiS, they persistently miss the intragap states. Addition- 
ally, correlations play an important role in NiO, and the real 
space approach is therefore more suitable to describe its com- 
plex electronic structure. Last but not least we wish not only 
to calculate the energy levels but the wavefunctions. 
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FIG. 1: Clusters used for the bulk (left) and the surface (right). Large 
dark spheres represent oxygen atoms, the central grey sphere the 
nickel atom, and the small spheres the surrounding effective core 
potentials. 



A. Doubly embedded cluster 



In order to model the system we use a doubly embedded 



cluster 



For the surface a NiOg^* cluster is used and for 
the bulk a NiOg cluster. Both clusters are first embedded 
in a shell of effective core potentials (ECP), that account for 
the electrostatic properties of the Ni ions in the immediate 
vicinity (Fig.[T]i. The whole structure is then again embedded 
in a charge point field (CPF) which describes the surrounding 
Madelung potential. 



B. Computational levels 

We use several correlation levels in the calculations, mainly 
for comparison and checking reasons, such as configuration 
interaction (CI) with single excitations (CIS), single and dou- 
ble excitations (CISD), single excitations with energy con- 
tributions from higher excitations [CIS(D)]^", quadratic con- 
figuration interaction (QCI), with singles (QCIS), doubles 
(QCISD), and doubles with higher excitation energy contribu- 
tions (QCISD-T). A method well suited to describe all (/-states 
is the multiconfigurational complete active space configura- 
tion interaction (MC-CAS or CASSCF), which can describe 
each state separately and can include dynamical correlations 
as well^' . TableUshows some results of NiOg^^ from different 
correlation levels. 



III. SPIN-ORBIT COUPLING 

The SOC, characterized by the interaction between the elec- 
tron spin moment and the orbital angular momentum, results 
in an additional effect of splitting of the electron energy levels. 
The high-quality linear and nonlinear (SHG) optical spectra of 
NiO given by Fiebig et alJ provide a good chance to directly 
compare the experimental and theoretical details of the fine 
structure. 



Using the CIS method and including SOC, we concentrate 
on the 21 lowest triplet states originating from the crystal field 
and spin-orbit splitting of the level of the Ni^"*" ion. Sym- 
metry analysis shows that SOC does not split the singlet states. 
It shifts, however, their energy due to the admixture of the 
triplet states. As a consequence some singlet-triplet transi- 
tions become allowed. 

Figure 2 in reference^ shows the calculated fine structure 
of bulk NiO. Including SOC leads to a slight lowering of the 
ground state energy. The splitting of the second excited state 
(Ai53p+= i5p4.— i?p+=71 meV) i'^Tt,) is in good agreement 
with the most accurate measurements (around 70 meV) of the 
optical absorption and SHG spectra of NiO^. In addition, our 



transition ' F 



3r+ 



^F^ at about 0.9 eV agrees well with the 



SHG structure around 1.0 eV by Fiebig et alJ . 

We have also found that the splitting of energies leads to 
additional electric dipole (ED) transitions, which corresponds 
to the SHG response that can be observed when SOC is taken 
into account. 



IV. SECOND HARMONIC GENERATION 
A. Theory 

When we consider a monochromatic plane wave interacting 
with a material, it induces a response which leads to a source 
term in the Maxwell equations. Expanding its inhomogeneous 
solution yields terms of the polarization (Ps), the magnetiza- 
tion (Ms), and the electric quadrupole (EQ) polarization {Qs), 
which are functions of the incident field^^-^'*. Taking into ac- 
count no more than one second-order transition and consid- 
ering the indistinguishability of the two incident photons, we 
have 
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In the right side, the two last upper indices denote the nature 
of the two absorptions, while the first one the nature of the 
emissions [e for ED, m for magnetic dipole (MD), and q for 
EQ induced transitions], and the lower indices in the subse- 
quent formulas represent the directions of the coiTesponding 
multipoles. To calculate x^'^^\ adopting the Coulomb gauge, 
neglecting the A'^ term (this term leads to a tensor completely 
analogous to the nonlinear optical conductivity tensor E*^ dis- 
cussed by Andersen et al?^) and setting n and a as the unit 
vectors in the propagation and the polarization direction of the 
light beam respectively, we get a perturbation term, a Taylor 
expansion of which keeping only the first two terms in i^n • r 
yields for the interaction elements: 
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where £ and B are magnitudes of the electric and magnetic 
fields, respectively, uj is the frequency of the incident light. 
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TABLE I: li-state and gap energies (in eV) at various correlation levels of the NiOg clustei^. 
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" Taken from Ref. 11. 



dfj, if/" and Q^j" are the components of the ED, MD, and 
EQ transitions along the directions indicated by the super- 
scripts. The microscopic formulas for x'^^'^K derived from 
the Liouville equation of motion, by applying second order 
perturbation theoryi^i2i2S, has already been given^^^ 
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where the indices /, j, k denote the nature of the transition 
(i.e., e for ED, m for MD, and q for EQ transitions) between 
states a and (3, and f{E) is the population distribution func- 
tion. The overbar means symmetrization with respect to the 
two incident photons. 



B. Symmetry analysis 



able to describe the translational invariance of the real crystal 
lattice and the different domains, which give different second 
harmonic signals 
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The NiOg^ cluster has a C4V symmetry with no inversion 
center, and therefore all rank three tensors are allowed, and 
their symmetrized forms are 
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Our cluster is in the ferromagnetic phase due to the model 
involving only one Ni atom. The real space approach is not 
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with three, one, and three independent matrix elements, re- 
spectively. Symmetrizing results in the vanishing of the xll™ 
element, in the Oh point group (bulk) we further have xVyz = 
"Xlz™ thus the whole tensor disappears. 

The tensors are of rank four when involving EQ transitions. 
Symmetrizing with respect to the last two indices can simplify 



the x'''^'^ tensor both for the bulk and for the surface, xfij^ki 
has nine nonvanishing contributions. For the group Oh we are 
treating with a T4 tensor^^ with three independent elements, 
the symmetrized form is 
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FIG. 2: Tensor elements of the bulk of NiO. Peak a corresponds to 



transition 1 A20 



with frequency lo, peak 6 to 1 A2g 



2 Tig with frequency lj, peak c to 1 A2g ^ 1' Tig with frequency 



peak d to 
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1 Tig with frequency 2lo. 



Furthermore we have also analyzed, using group theory, the 
contributions to the 'x^'^^') arising from magnetic domains ro- 
tated with respect to each other^^. In addition, a similar sym- 
metry analysis has been performed for ferromagnetic systems 
to obtain the nonlinear susceptibility tensor by Andersen et 
al. where even and odd contributions are attributed to spin- 
dependent transitions^^. 



C. Results 



Altogether there are three types of independent tensor el- 
ements, namely xlyT' Xxy{xyy ™d Xp^w^^^ in bulk NiO. 
The most distinct peak at 0.94 eV (Fig. |2| due to the MD 
transitions perfectly coincides with the experimental results^. 
Peak d, which arises from an EQ transition, also agrees very 
well. Further splitting of these peaks can easily be attributed 
to SOC and/or phonons. 

Considering certain experimental conditions, we repopulate 
the states in two ways^', that is, an energetically balanced way 
which consists of an equal distribution of electrons in all sub- 
levels with a given energy, and an energetically nonbalanced 
way where anisotropy is included at the same time, as demon- 
strated in Figures[3]and|4] The first case could happen when a 
system is prepared with an unpolarized pump laser pulse, and 
the results show that higher kinetic temperatures would result 
in excitations to different states. In the second way, which 
could happen with a polarized pump beam, x^x™ shows the 
same behavior as in the thermal distribution, while xlyT ex- 
hibits additional symmetry violating features. 
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FIG. 3: Energetically balanced repopulation of the NiO^ cluster 
and its effect on the Xij/T tensor element. The outer solid lines (black 
and red) represent the tensor element for fully populated ground and 
first excited state, the inner solid line (green) the suppressed peak 
when the population ratio equals 3:1. 
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FIG. 4: Energetically nonbalanced repopulation of the NiOjj"^" clus- 
ter and its effect on the xJyT and Xya:™ tensor elements. 



V. PHONONS 
A. Frozen-phonon approximation 

Up to now, we supposed a non-distorted lattice, which of 
course is not true in the real system. However, the system 
deviates from the full Oj, symmetry through the static lattice 
distorsion along the (111) axis due to spin ordering or dy- 
namic distortions like phonons. Here we consider the latter 
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FIG. 5: In the middle panel are the calculated levels for Oh symmetry, 
on the left for C4v and on the right for C2v The insets show the 
respective clusters (with the effective charge points). 
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FIG. 6: x'^"'^ and x"'""' tensors for the NiOg in the case of the F 
point of one optical phonon along the z direction (C4v symmetry). 
Note that peaks a and e are shifted with respect to each other by 5 
meV, a splitting which is visible only if zoomed in. 



case, which takes place both in the bulk and the surface. We 
treat the electronic movements quantum mechanically within 
the Born-Oppenheimer approximation while the motion of the 
nuclei is considered to happen on a much slower pace, giving 
time to the electrons to adiabatically follow them. 

The optical phonons lower the local symmetry of the 
NiOg^^" in our rocksalt structure, even at the T point. There 
are a few cases, however, where the local geometry on any Ni 
site is the same throughout the crystal, the level splittings of 
which are shown in Fig.|5] 

Now SHG becomes possible even within the ED approxi- 
mation since the inversion symmetry is broken. Fig.|6]depicts 
the x'^'^'^ and the 5(<=^™ tensor elements in the NiOg^^*' clus- 
ter for a phonon coordinate of 3% dislocation, at the T point, 
when the symmetry on every Ni site is C4V. Here, the energy 
splitting of the peak near 1 eV of the 3^'='='" (about 5 meV) 
is very close to the experimental splitting of the first intense 
SHG peak observed by Fiebig et alJ . 

A lower symmetry will lead to an even richer structure^l. 
Note that generally the triply degenerate states get split into 
three levels, but two of which are always almost degenerate 
in energy, and thus the peaks get split into two peaks and not 
three, as shown in Fig.|2] 
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FIG. 7: Splitting of the levels vs. displacement in the case of the 
transversal plane wave (see text) along the z direction (C2u symme- 
try). The splitting is an even function of q. 



B. Time evolution 

The calculations show that the relation between x^^^Hl) 
and q is not linear (Fig.|2ll. Nonetheless, since in our approx- 
imation we consider SHG to be instantaneous we calculate 
the time average over a whole period of (x'-^'^-')^ as shown 
in Fig. [8] We consider a simple harmonic model (uncoupled 
phonons), and take snapshots for amplitude of up to 3% lattice 
displacement with a step of 0.25%. It is found that the peaks 
do not get broadened but split, and the intensity of the peaks 



has an almost quadratic relationship to the amplitude of the 
displacement. The time resolved results shown in Fig. |9] and 
Fig. [To] exhibit the relationship between the tensors elements 
ixTxz ^iid xTyz) ^f"! the phononic coordinate q). Both the 
behaviors are clearly nonlinear, and one can distinguish three 
regimes^'. A similar behavior can be observed for the x™yz 
tensor element. 



FIG. 8: Time averaged | x'"='= | and {x"""" \ tensors for the NiOg ^° in 
the case of the F point of one optical phonon along the z direction 
(C4v symmetry) with an lattice displacement of amplitude of 3%. 
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FIG. 9: xSafz vs. photon energy and vs. q. A peak starts appearing 
after q = 0.25% at energy ~ 0.95 eV. 



VI. SPIN MANIPULATION 

A. Simple models for laser-induced dynamics 

We start our investigations of the spin manipulation with 
the simplest two-level and three-level systems^'. To obtain 
a spin mixed state, two excited states of definite (and oppo- 
site) spin are needed, which are mixed through SOC. Thus, 
the minimum scenario for optically induced spin dynamics is 
a four-level system: two almost degenerate pure spin states 
at low energy (ground states) and two also almost degenerate 
states at high energy which can be mixed by SOC-^. 



B. Ultrafast magnetic switching in NiO 

Based on the above suggested four-level system, we present 
an ultrafast magneto-optical switching mechanism in NiO. 
Fig. [TT| shows the population evolution of the lowest levels 



FIG. 10: xTyz vs. photon energy and vs. q. Two new peaks appear 
after q = 0.25%, one at energy ~ 1.590 eV which changes its phase 
at q = 2.0%, and one at energy ^ 3.65 eV. 




Time (fs) 

FIG. 11: Top: population evolution of the lowest lying levels in bulk 
of NiO. Bottom: Laser profile. 



in bulk NiO and the corresponding laser profile. In order to 
achieve switching we have to choose a suitable mixed-spin ex- 
cited state. We find that the maximum efficiency occurs when 
the matrix elements from both ground states have almost the 
same absolute values. The phase evolution of the different 
states governs the direction of the population transfer during 
the pulse. That is why Fig.fTTIexhibits a double-peak structure 
for the intermediate excited states. 

Fig.[T2]shows the dependence of the spin magnetic moment 
of the NiO(OOl) surface with respect to the duration and inten- 
sity of the pulse, after the pulse is switched off and the system 
has reached a static state. The static B-field is perpendicular 
to the surface and the laser pulse is linearly polarized at a di- 
rection parallel to the surface. Contrary to previous works^^, 
taking into advantage of the breaking of the time reversal sys- 
mmetry due to the B-field, we find optimal switching condi- 
tions for linearly polarized light, more specifically with prop- 
agation direction parallel to the static B-field — this is consis- 
tent with the experimental findings of Koopmans et al}^ in an 
experimental setup similar to ours. Our calculations indicate 
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FIG. 12: Contour plot of the final magnetic state of the NiOJ^ clus- 
ter after the application of a linearly polarized sec^^-shaped laser 
pulse and an external static magnetic field Bcxt — — 12.5 A/m 
vs. duration and intensity of the pulse. 

that the use of circularly polarized light can selectively acti- 
vate one channel of the A-process (see inset in Fig.fTTTi. only 
allowing either excitation or de-excitation, while a superposi- 
tion of the two polarizations can selectively control the per- 
centage of the population transfer^"^. Similar to our results, 
Stanciu et al. observed subpicosecond magnetization reversal 
in a ferrimagnetic system, however, with a circularly polarized 
laser pulse without magnetic field^^. 

Note that breaking of the time-reversal symmetry leads to 
a differentiation of the roles of (L-^) and {L). While the for- 
mer describes crystalline anisotropy and survives unharmed 
the existence of the external B-field, the latter one lifts spa- 
tial degeneracies and gives rise to dynamic effects of mag- 
netic origin (e.g. Kerr and Faraday effects). Thus linear and 
quadratic coupling to the (anti)-ferromagnetic order parame- 
ters become detectable in optics. 

C. Quantized phonons 

In order to consider phonons in a fully quantized picture 
we first calculate the normal modes of our cluster by diag- 
onalizing the dynamical matrix. Then we quantize the nor- 
mal modes and subsequently calculate the phonon-electron 

interaction by fully quantum-mechanically considering ^gy^, 
where q is the phononic coordinate as seen in Figures |9] 
and II Of ^. Then we include them together with the SOC and 
the external static magnetic field in one diagonalization step 
and therefore to eliminate the problem of applying twice a 
perturbative method on the wavefunction. Fig. [13] shows the 
magnetic state of the NiOg^^° cluster in the presence of one 
F-point optical phonon after the application of a linearly po- 
larized laser pulse and an external static magnetic field vs. du- 
ration and intensity of the pulse. The existence of optical 



FIG. 13: Contour plot of the magnetic state of the NiOg" cluster in 
the presence of one F-point optical phonon after the application of 
a linearly polarized sec^^-shaped laser pulse and an external static 
magnetic field Bext = — 12.5 A/m vs. duration and intensity of 
the pulse. 

phonons moves the switching parameters towards higher in- 
tensities and longer times, while decreasing the tolerance of 
the system for achieving switching. 

Temperature analysis indicates that the elevated lattice tem- 
perature makes the effect of switching more efficient, while 
electronic temperature (loosely defined) does not significantly 
alter the process, since the first one significantly changes sym- 
metry, while the second one does not^^. 



VII. SUMMARY AND OUTLOOK 

In summary, we have presented a quantum chemistry 
method for calculating different contributions to SHG both 
for the bulk and for the (001) surface of NiO, without involv- 
ing any empirical parameter. We consider not only the ED 
approximation but MD and EQ induced transitions as well, 
and we also calculate the phononic effects within the frozen- 
phonon approximation. It is shown that all of the following, 
the SOC, the surface, the nonlocal contributions of the laser 
field, and the transient lattice distortions, can give rise to peaks 
in the centrosymmetric bulk NiO. The results also indicate the 
possibility to disentangle the surface (001) from the bulk sig- 
nal due to different geometrical and frequency dependencies, 
which has been confirmed in a recent experiment'". 

We present a fully ab initio ultrafast magnetooptical switch- 
ing mechanism in NiO and show that controlled switching is 
possible by including SOC in order to take advantage of a A- 
process. The role of phonons is further discussed in a fully 
quantized picture as a symmetry-lowering mechanism in the 
switching scenario. 

The role of symmetry analysis is definitely significant 
through the whole paper. It has been used to classify the elec- 
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tronic states and judge the allowed dipole transitions in SOC, 
simplify the tensor elements of the high order transitions in 
SHG, and help us recognize the allowed transitions in the re- 
alization of ultrafast switching. 

Up to now, we only have considered one magnetic center 
(one Ni atom) in our investigations. Although the present 
single-site model can mostly account for AF ordering in NiO 
(both bulk and surface) because of the high spin density on the 
Ni sites, it is still necessary to involve more magnetic centers 
to study more subtle systems for possible applications such as 
magnetic logic units. 

Note that since our calculations were performed on a one- 
magnetic-center cluster we therefore did not need to consider 
gauge including atomic orbitals (GIAOs). However, clusters 
containing more than one active magnetic center should (per- 
haps) include the phase factor e*^ "" in the atomic orbitals. 

Although spin is not a good quantum number after the in- 



clusion of SOC the breaking of the time-reversal symmetry 
does not inhibit the calculation of its expectation value. How- 
ever, it is necessary to account for it in a full spinor basis (two- 
component) even though rotational invariance is violated in 
the cluster. For systems where the spin density is localized 
on one magnetic center only (i.e. NiO) one can consider this 
center as the origin of the spatial symmetry operations, but for 
multi-centered clusters the necessity of fully accounting for 
the spinors becomes imminent. 
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